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Nonlinear  Robust/ Adaptive  Controller  Design  for  an 
Air-breathing  Hypersonic  Vehicle  Model 
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This  paper  describes  the  design  of  a  nonlinear  controller  for  an  air-breathing  hypersonic 
vehicle.  To  overcome  the  analytical  intractability  of  this  model,  a  nominal  control-oriented 
model  is  used  to  derive  the  control  law.  The  nominal  model  has  unstable  zero  dynamics 
with  respect  to  the  output  to  be  controlled,  namely  flight  path  angle  and  velocity,  and 
presents  intricate  couplings  between  the  engine  dynamics  and  flight  dynamics.  The  flexible 
effects  have  not  been  included  in  the  analysis  yet.  Adaptive  control  techniques  and  robust 
control  techniques  are  applied  to  achieve  tracking  of  velocity  and  altitude  trajectories. 
Simulation  results  are  provided  to  show  that  the  derived  control  law  allows  to  achieves 
excellent  tracking  performance  on  the  nominal  model. 


I.  Introduction 

The  design  of  guidance  and  control  systems  for  air-breathing  hypersonic  vehicles  is  a  challenging  task,  due 
to  the  high  complexity  of  the  equations  of  motion,  stemming  form  strong  couplings  between  propulsive  and 
aerodynamic  effects,  significant  flexibility  associated  with  the  slender  geometries  required  for  these  aircraft, 
and  large  model  uncertainty.1-3  For  linearized  versions  of  hypersonic  vehicle  models,  several  results  are 
available  in  the  literature,1,2,4-8  which  consider  model  of  various  complexity.  Tournes  et  al.  employed  a 
nonlinear  variable  structure  control  approach  in  Ref.  9,  while  several  other  nonlinear  control  approaches 
have  been  proposed  for  non-flexible  models  of  more  general  types  of  vehicles.10  13  More  recently,  in  Ref.  14 
and  15,  linear  controllers  with  explicit  constraints  on  the  inputs  are  derived  for  linearized  models  using 
on-line  optimization  and  anti- windup  techniques. 

Especially  when  flexibility  effects  are  taken  into  account  in  the  model,  the  equations  of  motion  may  be¬ 
come  exceedingly  complex,  due  to  intricate  interactions  between  the  structural,  aerodynamic  and  propulsion 
systems.  As  this  happens,  models  derived  from  first  principles  can  be  used  effectively  only  for  simulations 
or  validation  purposes,  whereas  the  control  design  must  be  performed  on  the  basis  of  control-oriented  of 
reduced  complexity  which  approximates  the  behavior  of  a  first-principle  model  to  a  satisfactory  degree. 

In  Ref.  16,  a  control  oriented  model  is  derived  for  the  longitudinal  dynamical  model  developed  by  Bolen- 
der  and  Doman3  using  curve  fits  calculated  directly  from  the  forces  and  moments  included  in  the  truth 
model.  Feedback  linearization  techniques  are  then  applied  to  the  control-oriented  model  to  derive  a  nonlin¬ 
ear  controller.  Even  though  this  approach  allows  to  achieve  good  tracking  performance,  due  to  the  model 
complexity,  it  inevitably  leads  to  very  complicated  expressions  for  the  hight-order  Lie  derivatives  of  the 
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Table  1.  States,  control  inputs,  and  controlled  output  for  the  vehicle  model 


State 

V 

h 

a 

0 

7 

Q 

Vi 

Vi 


Input 


Vehicle  velocity 

Altitude  Ad 

Angle  of  attack  Se 

Pitch  angle  Sc 


Flight  path  angle,  7  =  0  —  a 
Pitch  rate 
i-th  Generalized  elastic  coordinate  V 

Time  derivative  of  r]i  h 


Fuel-to-air  ratio 

Diffuser  area  ratio 

Control  surface  deflection  (elevator) 

Control  surface  deflection  (canard) 


Vehicle  velocity 
Altitude 


Output 


regulated  outputs,  and  ultimately  for  the  controller.  Furthermore,  it  is  not  possible  to  perform  a  robustness 
analysis  when  model  uncertainty  is  taken  into  consideration. 

In  this  study,  we  take  a  first  step  toward  a  more  general  nonlinear  design,  not  based  on  dynamic  inversion, 
for  a  given  hypersonic  vehicle  model.  While  a  direct  application  of  nonlinear  control  system  design  techniques 
such  as  backstepping17  and  forwarding  18, 19  is  not  applicable  -  since  the  vector  field  structure  is  neither  in 
strict  feedback  form,  nor  in  feed-forward  form  -  some  basic  ideas  from  these  techniques  are  used  to  derive  the 
controller.  Following  a  classic  sequential  loop  closure  approach,20  the  model  dynamics  have  been  decomposed 
into  subsystems  each  of  which  has  its  own  control  input.  In  particular,  the  thrust  is  used  to  control  the 
velocity  of  the  vehicle,  the  angle  of  attack  is  used  to  control  the  flight  path  angle  and  the  altitude,  and  the 
moment  is  then  used  to  control  the  angle  of  attack/pitch  rate  subsystem.  In  each  subsystem  the  control 
law  is  derived  using  a  “mixed”  approach  based  on  adaptive  techniques21  and  robust  techniques  such  as 
the  small-gain  theorem  for  input-to-state  stable  systems.22,23  Since  dynamic  inversion  is  not  applied,  the 
controller  does  not  depend  explicitly  on  the  model  coefficients,  and  can  in  principle  be  designed  to  provide 
robustness  with  respect  to  parameter  uncertainties  and  variations  within  a  given  range. 

The  vehicle  model  considered  in  this  work  is  the  assumed-modes  version24  of  the  model  of  the  vehicle 
longitudinal  dynamics  developed  By  Bolender  and  Doman.3,25  Because  of  the  complexity  of  this  model, 
a  control-oriented  model  has  been  developed  following  the  approach  used  in  Ref.  16.  The  outputs  to  be 
controlled  are  the  vehicle  velocity  and  altitude.  The  control  inputs  are  the  elevator  and  canard  deflections, 
and  the  fuel-to  air-ratio  of  the  scramjet  engine.  For  this  preliminary  result,  the  presence  of  the  flexible 
dynamics  is  not  taken  into  account  directly  at  the  design  level,  but  regarded  as  a  perturbation  on  a  nominal 
model,  and  its  effect  evaluated  in  closed-loop  simulations. 

The  paper  is  organized  as  follows:  in  Section  II  the  first-principle  model  is  briefly  introduced,  and 
simplified  curve- fitted  and  control-oriented  models  are  presented;  the  control  law  design  is  described  in 
Section  III,  and  simulation  results  for  the  curve- fitted  model  are  given  in  section  IV.  Finally,  Section  V 
concludes  the  paper  with  a  brief  summary  of  the  results,  and  a  discussion  of  further  work. 

II.  Vehicle  model 

The  vehicle  model  considered  in  this  paper  is  the  assumed-modes  version24  of  the  model  originally  devel¬ 
oped  by  Bolender  and  Doman3, 25  for  the  longitudinal  dynamics  of  a  flexible  air-breathing  hypersonic  vehicle. 
The  equations  of  motions,  derived  using  Lagrange’s  equations  and  compressible  flow  theory,  include  flexi¬ 
bility  effects  by  modeling  the  vehicle  as  a  single  flexible  structure  with  mass-normalized  mode  shapes.  The 
scramjet  engine  model  is  taken  from  Chavez  and  Schmidt.26  Since  this  first-principle  model  is  analytically 
intractable,  a  control-oriented  version  of  the  model  has  been  derived  by  replacing  complex  expressions  of 
the  aerodynamic  forces  and  moments  with  curve-fitted  approximations,16  and  by  retaining  the  first  three 
flexible  modes.  In  this  preliminary  work,  the  control-oriented  model  is  used  in  this  paper  for  both  control 
design  and  simulation. 
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A.  Assumed-Modes  Model 


The  assumed-modes  model24  considered  in  this  paper  yields  the  equations  of  motion  of  the  vehicle  longitu¬ 
dinal  dynamics  written  in  the  stability  axes  as 


V 

h 

a 

0 

Q 


Vi 


T  cos(o)  —  D 
m 

V  sin (0  —  a) 

L  +  T  sin(o') 
mV 


g  sin(#  —  a) 


+  Q  +  ^  cos (0 


Q 


a) 


M 

lyy 

2 Ci^iVi  Vi  T  ,  i  —  1,  2,  3 


(1) 


where  m  is  the  vehicle  mass,  Iyy  the  moment  of  inertia,  g  the  acceleration  of  gravity,  and  Q  and  uji  are  the 
damping  factor  and  the  natural  frequency  of  the  z-th  flexible  mode,  respectively.  The  aerodynamic  forces 
and  moments  enter  (1)  as  lift  L,  drag  D,  thrust  T,  pitching  moment  M  about  the  body  y-axis,  and  the 
generalized  elastic  forces  A^,  i  =  1,2,3.  The  model  comprises  five  rigid-body  state  variables  xr,  six  flexible 
states  77,  and  four  control  inputs  u,  as  defined  in  Table  1: 

X  =  [v,  a,  Q,  h,  e]  ,  T)  =  [j?i,?7i,  772,772, 773,773]  ,  u=[$,Ad,Se,5c\ 

whereas  the  output  to  be  controlled  is  selected  as  y  =  [V,  h] ' .  The  equations  that  define  the  aerodynamic 
and  generalized  forces  and  moments  (not  reported  here  for  reasons  of  space  limitation)  are  quite  intricate 
functions  of  the  state  variables,  the  input  and  the  plant  parameters,  which  do  not  admit  a  closed-form 
representation.  The  interested  reader  is  referred  to  the  cited  references  for  details. 


B.  Curve-Fitted  Model  and  Control-Oriented  Model 

The  first-principle  model  (1)  leads  to  equations  of  motion  in  which  the  relationships  between  control  inputs 
and  controlled  outputs  are  not  explicit.  For  controller  design,  a  simplified  model  has  been  developed  fol¬ 
lowing  the  approach  used  in  Ref.  16.  The  simplified  model,  referred  to  as  the  curve- fitted  model  (CFM), 
approximates  the  behavior  of  the  assumed-modes  model  by  replacing  the  aerodynamic  and  generalized  forces 
and  moments  with  curve- fitted  approximations.  The  resulting  non-linear  model  offers  the  advantage  of  being 
analytically  tractable  (albeit  still  quite  complex)  and  more  suitable  for  control  design,  while  retaining  the 
relevant  dynamical  features  of  the  first-principle  model.  The  approximations  of  the  forces  and  moments 
employed  in  the  CFM  are  given  as  follows 

T  «  +  CAd,<s>(®)Ad  +  Cr,o(oO  +  CT^V 

M  «  ZTT  +  <pSCM{cL,5e,8c,rj) 

L  «  qSCL(a,  Se,  Sc,  rj) 

D  «  qSCD(a,  Se,  Sc,  rf) 

Ni  «  Nf  a2  +  N?a  +  N?'5e  +  N?‘6C  +  JV°  +  NVV 
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where  the  thrust-to-moment  coupling  coefficient  zt  and  the  mean  aerodynamic  chord  c  are  given  constants, 
q  denotes  dynamic  pressure,  and 


CT,Ad(a) 

Ct,o 

Cm{ol ,  Se,  5C,  rf) 
CL(a,5e,5c,r]) 


CD{a,Se,Sc,r]) 


C|"3  a3  +  C£aV  +  C£“a  +  C$ 

C^d0?  c?  +  Cyd“2  a2  +  C^d<xa  +  C^d 
cf  a3  +  Cfa2  +  C%a  +  <7 
Cm  a2  +  Cm  a  +  C^Se  +  C^Sc  +  C%  +  C^ri 
C£a  +  Ci‘6e  +  Ci‘5c  +  C°L  +  C> 

Cf)  a2  +  CpCt  +  C q  52  +  Cp  8e  +  <52  +  Cp  5C  +  £?£>  +  C^fy . 


In  this  study,  the  CFM  is  used  primarily  for  simulation,  while  for  the  development  of  the  controller  a 
control-oriented  model  (COM)  is  employed,  obtained  from  the  CFM  by  neglecting  the  flexible  dynamics.  As 
a  result,  the  flexible  dynamics  are  not  taken  into  account  directly  at  the  design  level,  but  are  considered  as 
perturbations  on  the  COM,  and  their  effects  evaluated  in  simulation.  Consequently,  the  COM  comprises  of 
the  five  rigid-body  state  variables  and  four  control  inputs.  Among  the  available  control  inputs,  the  canard 
surface  deflection  (5C)  will  be  used  to  decouple  the  lift  force  from  the  elevator,  the  fuel  to  air  ratio  (<F)  will 
be  used  to  command  thrust,  and  hence  to  control  the  vehicle  velocity,  while  the  elevator  surface  deflection 
will  ultimately  be  used  to  control  the  angular  dynamics  through  the  pitch  moment.  Since  the  diffuser  area 
ratio  (Ad)  will  not  be  used  as  independent  control  input,  its  value  will  be  set  to  Ad  =  1  for  the  remainder  of 
the  paper.  Furthermore,  to  facilitate  the  development  of  the  controller,  the  flight-path  angle  7  =  0  —  a  will 
be  used  in  place  of  the  pitch  angle  as  a  state  variable.  The  equations  of  motions  for  the  COM  are  given  by 


where 


V 

h 

7 

a 

Q 


T  cos  a  —  D 

- asm  7 

m 

V  sin  7 


L  +  T  sin  a 
mV 

L  +  T  sin  a 
mV 


9 

-cos  7 
+  Q  +  77  COS7 


M 

lyy 


T  ph  Ct, &(&)*&  +  Ct,  o(^) 

M  pz  ztT  +  qcSCM(&,  <5C,  0) 
L  «  qSCL(a:5e:5c,  0) 

D  «  qSCD(a,6e,6c,0) 


(2) 


(3) 


and  the  expression  of  Ct$(ol)  reads  as 


CT, 0  =  (C^da  +  C^)a3  +  (C^da  +  C^)a2  +  (C$d%  C$)a  +  C$d  #  C %  . 


III.  Robust / Adaptive  Controller  Design 

The  control  problem  considered  in  this  paper  is  to  track  given  velocity  and  altitude  references,  generated 
using  filtered  step  signals  to  account  for  physical  limits  of  performance  of  the  vehicle.  The  approach  we 
consider  for  control  design  aims  at  avoiding  writing  the  model  in  the  so-called  normal  form  for  the  purpose 
of  inverting  its  dynamics,  as  traditionally  done  in  feedback  linearization  schemes.  The  first  reason  behind 
this  choice  lies  in  the  fact  that  the  model  is  non-minimum  phase  with  respect  to  the  input  [3>,5e]  and 
the  output  to  be  controlled  [V,  h\  .  Apart  from  a  weakly  non-minimum  phase  behavior  stemming  from  the 
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flexible  effects  for  the  CFM,  both  the  CFM  and  the  COM  possess  an  exponentially  unstable  zero-dynamics 
(with  respect  to  the  given  I/O  pair)  due  to  the  coupling  between  elevator  deflection  and  lift.  The  occurrence 
of  an  unstable  zero-dynamics  renders  a  design  based  on  dynamic  inversion  problematic,  and  may  lead  to 
poor  closed-loop  performance  or  even  instability.16  Secondarily,  even  a  normal  form  based  on  the  COM 
leads  to  a  very  complicated  expression  for  the  high-order  Lie  derivatives,  which,  together  with  the  fact  that 
a  closed-form  expression  of  the  diffeomorphism  is  not  available,  render  a  robustness  analysis  intractable.  For 
these  reasons,  we  propose  in  this  work  a  nonlinear  design  based  on  classic  sequential  loop-closure  arguments, 
which  uses  a  combination  of  adaptive  and  robust  techniques  and  avoids  resorting  to  dynamic  inversion.  The 
methodology  proceeds  by  decomposing  the  equations  of  motions  of  the  COM  into  functional  subsystems, 
which  are  controlled  separately  using  suitable  virtual  control  inputs.  For  each  subsystems,  an  appropriate 
control  law  is  designed  by  taking  into  consideration  the  structure  of  the  equations,  the  level  of  interaction 
with  the  rest  of  the  dynamics,  and  the  time  scale  of  the  variables  to  be  controlled.  For  the  slower  portion 
of  the  dynamics  (velocity  and  altitude),  first  a  control  law  with  adaptive  drag  compensation  is  derived  for 
the  velocity  subsystem  by  controlling  the  thrust  from  <f>.  Then,  a  controller  for  the  altitude  dynamics  is 
developed  using  the  flight-path  angle  as  a  virtual  reference  input.  Finally,  the  fast  angular  dynamics  with 
states  [a,  Q\  is  used  as  a  servo-controller  (commanded  by  the  pitch  moment)  to  let  j(t)  track  the  required 
virtual  reference  command. 

Since  the  COM  is  obtained  from  a  curve-fitted  approximation  of  a  first-principle  model,  it  is  fundamental 
that  the  control  law  provides  robustness  with  respect  to  uncertainty  on  the  parameters  of  the  plant  model. 
In  developing  the  controller,  it  is  assumed  that  all  the  coefficients  of  the  COM  (apart  from  obvious  ones 
corresponding  to  physically  measurable  quantities  or  known  constants)  are  subject  to  a  possibly  large  un¬ 
certainty.  In  particular,  the  design  must  ensure  that  the  control  objectives  (in  terms  of  boundedness  of  all 
internal  variables  in  closed-loop  and  convergence  of  the  tracking  error)  are  maintained  for  all  possible  values 
of  the  plant  parameters  and  for  all  initial  conditions  ranging  over  given  compact  sets. 


A.  Adaptive  Controller  for  the  Velocity  Subsystem 

Substituting  the  expression  of  the  thrust  T  in  Eq.  (3)  into  the  first  equation  of  the  COM  dynamics  (2),  the 
equation 

mV  =  Ct,$  (o)  cos  a  <f>  +  Ct, o  (ck)  cos  a  —  D  —  mg  sin  7 

is  obtained  for  the  dynamics  of  the  vehicle  velocity,  in  which  the  input  <f>  appears  explicitly.  Since,  due  to  fuel 
consumption,  the  mass  of  the  vehicle  is  uncertain  and  changes  with  time,  we  set  m  =  m(t)  •  mo,  where  mo 
is  the  nominal  mass  value  and  m(t)  is  a  multiplicative  uncertainty.  Since  fh(t)  changes  slowly  with  respect 
to  the  time  scale  of  the  references  to  be  tracked,  it  is  considered  constant  during  each  setpoint  tracking 
maneuver.  Defining  V  as  the  difference  between  V  and  the  reference  Vref ,  the  dynamics  of  the  tracking  error 
for  the  vehicle  velocity  is  written  as 


mV  =  Ct,<p(cv)  cos<a  <F  +  Ct, o(ck)  cos<a  —  D  —  mg  sin  7  —  mVre f. 

By  introducing  the  vector  of  uncertain  parameter 

n  _  r /o<E>od  ,  p-y a.3  siAdCt2 ,  ^ia2  s^AdOn,  p-y A&\  riO 

yJ\  —  L°T  5  ,  O^p  ,  Kjrp  ,  L yji  “r  ,  w^p  ~c  ?  O^p  ~c  -p  Op->, 

SCg2 ,  S'Cg ,  SC§ ,  SC% ,  SC g ,  SC5J ,  SC°D ,  rh] ' 

and  defining  the  regressor  and  the  input  matrix 

(£,?/)  =  [0, 0,  0,  0,  —  a3  coscq  —  a2  coscq  —  <acoscq  —  coscq 

qa2,qa,q5l,q5e,q6l,q5c,q,  m0g  sin  7  -  m0Vret } 

B1(a)  =  [a3  cos  a,  a2  cos  a,  a  cos  a,  cos  a,  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0] 

equation  (4)  can  be  written  in  the  linearly  parameterized  form 


(4) 


mV  =  01B1(a)^  —  Tx(x, u)0\ .  (5) 

Functional  controllability  of  the  vehicle16  implies  that  01Bi(a)  ^  0  for  all  values  of  a  within  the  feasible 
envelope  of  flight  conditions.  This  conditions  is  assumed  to  hold  for  all  0 1  within  a  given  compact  set 
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@1  C  M16  chosen  to  represent  all  possible  variations  of  the  considered  plant  parameters.  The  next  step  is 
to  derive  a  certainty-equivalence  controller  using  the  Lyapunov  function  candidate  W\(V)  =  ^V2 .  The  Lie 
derivative  of  W\(V)  along  the  vector  field  defined  by  the  right-hand-side  of  Eq.  (5)  is  given  by 

Wi(V)  =  V  (^Si(a)$  -  ®i(®,w)0i)  . 

Choosing 

where  k\  >  0  is  a  gain  parameter,  the  derivative  of  the  Lyapunov  function  candidate  is  rendered  negative 
definite.  Replacing  the  uncertain  parameter  with  a  suitable  estimate  vector  #i,  the  input  <f>  is  ultimately 
selected  as 

®  =  I^—A-kiV  +  <S>'1(x,u)e1).  (6) 

OMa) 

Finally,  letting  6\  =  6\  —  #i,  the  Lyapunov  function  derivative  is  written  as 

m  =  V  (§[  -  01)B1(a)^  -  V  ®i(®,  w)(0i  -  0\) 

=  -k1V2  +  e'1(V^1(x,u)-VB1(a)^)  (7) 

The  choice  of  the  update  law  for  the  parameter  estimate  #1,  aiming  at  canceling  the  sign-indefinite  terms 
in  Eq.  (7),  will  be  postponed  till  later  in  this  section.  Note  that  to  ensure  non-singularity  of  the  control 
law  (6)  over  the  envelope  of  feasible  flight  conditions,  it  suffices  to  constrain  the  estimates  0\(t)  to  evolve 
within  a  suitable  compact  convex  set  which  includes  the  compact  set  ©i.  This  can  be  easily  accomplished 
by  resorting  to  smooth  parameter  projection  methods  in  the  definition  of  the  update  law.21 

B.  Robust  Integral  Controller  for  the  Altitude  Subsystem 

Given  a  desired  reference  profile  hre f  for  the  altitude,  the  corresponding  tracking  error  is  defined  as  h  = 

h  —  href.  Since  for  small  values  of  the  flight  path  angle,  h  =  Esin 7  «  Eq ,  it  follows  that 

fl  —  fl  href  ~  7  ^  href  . 

Choosing  a  reference  command  for  the  flight-path  angle  as 


Aref  — 


(h  href)  T  href 
V 


(8) 


where  kh  >  0  is  a  gain  parameter,  the  corresponding  dynamics  for  the  altitude  tracking  error  satisfies 


h  =  —khh 


and  thus,  if  the  flight-path  angle  is  controlled  to  follow  the  command  7ref,  the  altitude  tracking  error  is 
regulated  to  zero  exponentially  fast.  Note  that  the  reference  trajectory  (8)  is  well-defined,  as  V  >  0  in  any 
feasible  flight  conditions.  Next,  we  proceed  to  design  a  control  law  to  let  y(t)  track  7ref(t)  asymptotically, 
using  the  angle  of  attack  as  a  virtual  control  input. 

In  controlling  the  flight-path  angle  dynamics  of  the  vehicle,  one  of  the  major  obstacle  is  its  non-minimum 
phase  characteristic,  due  to  the  explicit  dependence  of  the  lift  force  on  the  control  surfaces.27  In  order  to 
eliminate  this  non-minimum  phase  behavior,  the  canard  is  used  to  suppress  the  influence  of  the  elevator  on 
the  lift  force  by  ganging  the  canard  deflection  to  that  of  the  elevator.  Specifically,  the  choice 


£c  = 


(9) 


achieves  direct  cancelation  of  the  elevator- to- lift  contribution.  To  account  for  parameter  uncertainty,  the 
selection  of  the  coupling  coefficient  between  the  canard  and  the  elevator  deflection  can  be  endowed  with  a 
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self-tuning  mechanism.  In  this  preliminary  result,  for  the  sake  of  simplicity,  the  choice  of  Sc  will  be  kept 
fixed  as  in  (9).  After  compensation,  the  lift  equation  reads  as 

L  =  qS(CZa  +  C°L).  (10) 

Using  Eq.  (10)  and  the  first  equation  in  (3),  the  dynamics  of  the  flight-path  angle  can  be  written  in  the  form 

mV 7  =  qS  (C^a  +  C£)  +  sin  a  +  Ct,o(&)  sin  a  —  mg  cosy  (11) 

Defining  the  tracking  error  for  the  flight-path  angle  as  7  =  7  —  7ref  ,  it  follows  that 

7  =  — ^  (qS  (Cl a  +  C£)  +  Ct&{ol)  sum  4>  +  Ct${ol)  shm  —  mg  cos(7  +  7ref)  —  mU7ref) 

•=  /O,  7,4,$,7ref).  (12) 


The  velocity  and  altitude  reference  trajectories  bring  the  vehicle  to  a  desired  trim  condition,  defined  as 

x*  =  [U*,h*,a*,7*,Q*]7,  u *  =  [$*,£*,£*]' 

where,  in  particular, 

U*  =  lim  Uref(t),  h*  =  lim  href(£),  7*  =  lim  7ref(t)  =  0,  Q*  =  0 

t — >00  t — >00  t — >00 


and  the  trim  condition  for  <1>  is  determined  by  the  steady-state  value  of  the  parameter  estimates  6*  = 
lim^oo  Qi  (t) 

$*  =  - - - - *',  (»*,«*)  ft-  (13) 

91  B\(a*) 

Furthermore,  we  let  q  and  a*  denote  respectively  the  value  of  the  dynamic  pressure  and  the  angle  of  attack 
at  the  desired  trim  condition.  Expanding  Eq.  (12)  around  ck*,  yields 

7  =  /o(a*,7,<7,$,7ref)  +  fi(a*,$)(a  -  a*)  +  R{a  -  a*)  (14) 


where 


A  (a*,*) 


1  r  not  i  d(Cr,s  (a)  sin  a)  d(CT, o(a)sina) 

m V[qS°L+  fa  |a=a**  +  ^ 


and  R(a  —  a*)  is  a  function  vanishing  at  a  =  a*  together  with  its  first  derivative.  To  obtain  an  equilibrium 
at  trim,  the  following  condition  needs  to  be  satisfied 


/o(a*,0,g*,$*,0)  =  0.  (15) 

Solving  (15)  yields  the  trim  condition  a*  that  enforces  the  equilibrium  in  7  =  0.  In  particular,  it  is  seen  that 

/o(a*,7,g,$,7ref)  =  (q-q*)S(CZa* +C°L)  +  CTAa*)s™a*($-$*)-mVjie{ 

- mg  cos(7ref)  +  mg  [(1  -  cos7)  cos(7ref)  +  sin7  sin(7ref)] 

which  vanishes  at  trim.  The  value  at  trim  for  the  angle  of  attack  must  necessarily  be  generated  by  the 
controller  to  maintain  the  desired  trim  condition  in  closed- loop.  However,  it  is  not  possible  to  compute  a* 
analytically  using  (15),  since  /o  depends  on  uncertain  plant  parameters.  To  overcome  this  problem,  the 
system  is  augmented  with  an  integrator,  whose  state  has  the  purpose  of  reconstructing  asymptotically  a*  in 
closed-loop.  The  integrator  dynamics  is  given  as 

i  =  v  (16) 

where  the  input  v  =  v(j)  is  chosen  as  a  function  of  7  which  vanishes  at  7  =  0.  Letting  £  =  £  —  a*,  the 
interconnection  of  system  (14)  with  system  (16)  reads  as 

7  =  /1  («*,$)(«  -  a*)  +  f0(a*,y,q,$,yle{)  +  R(a  -  a*) 

I  =  v(j)  (17) 
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where  the  angle  of  attack  (which  plays  the  role  of  a  virtual  control  input  for  system  (17))  and  the  function 
^(7)  must  be  chosen  in  such  a  way  that  the  equilibrium  [7,  £]  =  [0,  0]  is  asymptotically  stable.  The  desired 
angle  of  attack  is  chosen  as 

ad=£-k  27 ,  (18) 

where  &2  >  0  is  a  gain  parameter.  Assuming  £(£)  —>  a*,  the  first  term  in  the  right-hand  side  of  equation  (18) 
enforces  the  equilibrium  at  trim,  whereas  the  second  term  has  a  stabilizing  effect,  since  /i(<a*,  3>)  0  in  the 

range  of  plant  parameters  and  state  variables  of  interest.  Therefore,  imposing  (18)  yields 

a  —  a*  =  £  —  &27  +  ol  ,  (19) 


where  a  denotes  the  tracking  error  for  the  desired  dynamic  of  the  angle  of  attack.  The  control  v  (7)  is  chosen 
simply  as 

^(7)  =  -/c27- 

To  proceed  with  the  analysis,  we  set  a  =  0  and  verify  if  there  exists  a  choice  of  the  gain  k<±  for  which  the 
system 


4  =  /o(V,7,9,$,7ref)+/l(V,$)(C>fc27)+i?(|-fc27) 

l  =  ^(7) 

is  asymptotically  stable.  To  this  purpose,  the  following  change  of  coordinates 

X  ^  -P  (  *  7 

/l  («*,$) 


(20) 


(21) 


is  applied  to  system  (20).  Since  the  velocity  dynamics  are  slower  than  the  angular  dynamics,  <h(t)  can  be 
considered  constant  in  the  time  scale  of  system  (20).  As  a  result,  it  is  possible  to  consider  the  following 
approximation  for  the  dynamics  of  the  integral  error  in  the  new  coordinates  y: 


X  ~  -*27  - 


1  i 

A(a*,$)7' 


For  the  sake  of  notational  convenience,  in  what  follows  we  suppress  the  arguments  from  the  functions 
/o(cF,  7,  g,  7ref)  and  /i(a7,<f>).  In  the  new  coordinates,  the  system  (20)  assumes  the  form 


X 

7 


-x  -  k2j  +  (k2  -  7)  7  -  7  _  jrR(x  -  (k2  -  J-)  7) 


fa  ~  (/i*2  -  1)  7  +  fiX  +  R(X  ~  (k2  -  — )  7)  • 


(22) 


Looking  at  the  linear  approximation  of  the  system  (22)  around  the  origin,  which  reads  as 


X 


' x  +  [~h 


—  mg  sin(7ref)]  7 
/ 1 


7  =  -[/i*2  -  1  ~mg  sin(7ref)]  7  +  fix , 


(23) 


it  becomes  clear  that,  when  7  =  0,  the  y-subsystem  in  Eq.23  is  exponentially  stable,  and  when  X  =  0 
the  7-subsystem  is  exponentially  stable  for  large  enough  (recall  that  /i(a7,<F)  ^>0).  As  a  result,  each 
subsystem  is  input-to-state  stable22  (ISS),  with  linear  gains  given  by 


712 


721 


sup 

t>0 


—  mg  sin(7ref(t)) 
/ 1 


sup 

t>0 


_ h _ 

fik2  -1  -  mg  sin(7ref(i)) 


The  gain  of  the  overall  interconnection  is  given  by  the  product  of  the  single  gains,  from  which  it  is  evident 
that  there  exists  a  value  >  0  such  that  712721  <  1  for  all  k^  >  k%-  As  a  consequence,  by  the  small-gain 
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theorem,  the  system  (23)  is  globally  exponentially  stable  for  V&2  >  k^-  As  a  consequence,  the  system  (22)  is 
locally  exponentially  stable  as  well.  It  should  be  noted  that,  while  the  linear  system  (23)  possesses  an  infinite 
gain  margin  with  respect  to  &2,  the  original  nonlinear  system  (22)  need  not  share  the  same  property,  due  to 
the  appearance  of  &2  in  the  higher-order  term  R(£  —  £27).  As  a  consequence,  the  choice  of  the  gain  k2  must 
be  made  to  achieve  the  largest  possible  domain  of  attraction  for  the  origin  of  system  (22),  for  all  possible 
values  of  the  plant  parameters.  This  can  be  accomplished  using  numerical  optimization  methods  by  ensuring 
that  that  the  Lie  derivative  of  the  Lyapunov  function  candidate  W2O7,  x)  =  \  (72  +  X2)?  evaluated  along  the 
vector  field  of  system  (22),  is  negative  definite  in  a  level  set  Qc  =  {VL2(7,  x)  <  c}  which  contains  the  required 
envelope  on  initial  conditions,  for  all  the  considered  range  of  the  values  of  the  plant  parameters.  Specifically, 
letting  z  =  [7,  x]  f°r  notational  convenience,  and  substituting  (19)  in  equation  (14),  the  (7,  x)-dynamics  can 
be  written  in  the  form 

z  =  <Pi(z)  +  <fi  (z,a)a.  (24) 

It  is  assumed  that  k 2  has  been  selected  such  that  the  system  i  =  ipi(z)  has  a  locally  exponentially  stable 
equilibrium  at  2  =  0,  and  satisfies  the  uniform  Lyapunov  stability  property23  with  respect  to  the  Lyapunov 
function  W2(z)  =  |||^||2.  In  particular,  it  is  assumed  that  there  exists  a  constant  c\  >  1  such  that  the  set 
QCl  =  {z  :  W2(z)  <  ci  +  1}  is  compact,  contains  all  the  initial  conditions  of  interest  and  is  such  that,  for  all 
values  of  the  uncertain  plant  parameters, 

dWo 

V  <  -T i(z)  Wzenci  (25) 

oz 

where  Yi  (z)  is  a  continuous  and  positive  definite  function.  Following  the  approach  considered  in  Ref.  23, 
the  Lyapunov  function  W2(z)  is  modified  as 


W2(z)  =  c% 


W2(z ) 

Cl  7  1  -  w2(z) 


so  that  W2(z)  is  proper  on  the  set  QCl  while  at  the  same  time  satisfying 


d  Wo 

<pi(z,0)  <  -Ti(z)  Vz  e  QC1 

for  some  continuous  and  positive  definite  function  Yi(^).  Moreover,  since  the  stability  is  exponential,  there 
exists  a  >  0  and  r  >  0  such  that  Yi (z)  >  ^||^||2  for  all  2:  such  that  \\z\\  <  r. 


C.  Robust  Controller  for  the  (a,  Q)-Subsystem 


Next,  we  consider  the  problem  of  controlling  the  angular  dynamics  to  regulate  asymptotically  a(t)  to  zero. 
For  this  purpose,  the  control  input  for  the  angular  dynamics  given  by  the  (a,  Q)-subsystem  in  Eq.  2  is 
initially  selected  as  the  pitch  moment  M.  The  approach  we  take  for  controller  design  is  based  on  robust 
semi-global  stabilization  methods.19,23 

Let  ks  >  0  be  a  gain  parameter,  and  consider  the  change  of  coordinates 


yielding  a  system  of  the  form 
where 


Tn  = 


[a,  Q]'-C=[Ci,C2]'  :=[«,  £Q]' 

(  =  k3F0(  +  G1M(,z)  +  G2M, 


0  1  ,  G1=  1  ,  G2  = 

00/  0  / 


(26) 

(27) 


lyyks 


and  cp2  ((,z)  is  a  function  vanishing  at  the  origin.  The  desired  value  for  the  pitching  moment  is  selected  as 

Md  =  —kl  a0  a  —  k3  a\  Q  (28) 

=  kl(  — a0  ( 1  —  ai  C2)  (29) 
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where  ao  and  a\  are  such  that  p( A)  =  A2  +  aiA  +  ao  is  a  Hurwitz  polynomial,  the  system  (27)  can  be  rewritten 
in  the  following  form: 

c  =  hFC  +  Gnp2(t,  z)  +  G2(M  -  Md)  (30) 

where 

°  i  \ 

F=  _ao_  _ai_ 

\  lyy  lyy  ) 

Since  by  construction  F  is  a  Hurwitz  matrix  for  all  Iyy  >  0,  for  any  positive  numbers  I\  <  I2  there  exists 
a  fixed  symmetric  and  positive  definite  matrix  P  that  satisfies  the  Lyapunov  equation  F  P  +  PF  <  —I  for 
all  Iyy  G  [ii,  I2].  Applying  the  change  of  coordinates  (26)  to  the  equation  (24),  it  is  possible  to  group  the  z- 
and  ("-system  equations  as 


i  =  <Pl(2)  +  £l(^Cl)Cl  (31) 

C  =  k3FC  +  G1<p2(C,z)  +  G2(M-Md)  (32) 

and  consider  the  overall  Lyapunov  function  W3(zX)  =  W2(z)  +  £  P(.  Since  W2  is  positive  definite, 
FjFz  \z=Q  =  0  and  therefore  FjPz  =  v2(z)  z  for  some  continuous  function  v2(z).  Moreover,  since  <p2(0,  0)  =  0, 
it  is  possible  to  write 

^2  ((,Z)  =  <p2il(C,z)C  +  <P2,2(C>z)z 

for  some  continuous  functions  p2p ((,z)  and  cp2}2(C,z).  Due  to  the  fact  that  the  function  W3(zX)  is  proper 
on  QCl  x  M2,  there  exist  constants  c2  >  0  and  C3  >  0  such  that  the  set  13  =  {(2,  £)  :  c2  <  W3(z,  £)  <  C3}  is 
compact,  and  the  level  sets  DC2  =  {(2,  £)  :  W3(zX)  <  c2}  and  DCs  =  {(£,(")  :  W3(zX)  <  C3}  are  such  that 
ClC2  C  6r,  and  DCs  contains  all  initial  conditions  of  interest.  Defining 

Mi  :=  max  \\v2{z)  Ci)|| 

(z,c)eu 

M2  :=  2||P||  ||Gi ||  max  11^2,1  (C,2) II 

0,C)eo 

Ms  :=  2||P||  ||G2||  max  ||9?2)2(C,^)II 

GX)eU 


it  follows  that  the  derivative  of  W3(z,  £)  along  the  trajectories  of  the  system  (31)  satisfies 


Ws  =  ^Mz)F^MzX1)CiP^[k3FCPG1p2(Cz)}P^G2(M-Md) 

o  dWo 

<  -'Ti(z)  -  h IICII2  +  -gf  Mz,Ci)Ci  +  2C  PGW2&Z)  +  2C  PG2(M  -  Md) 

<  -Ti(z)  -  fc3|C|2  +  Mi \z\ ICI  +  M2ICI2  +  MskllCl  +  2C' PG2(M  -  Md) 


for  all  (z,£)  inQC3,  and  for  all  values  of  the  plant  parameters  in  a  given  compact  set.  As  a  result,  there 
exists  &3  >  0  such  that  for  all  k%  >  k%  the  trajectories  of  system  (31)  when  M  =  M d  are  captured  by  the  set 
DC2.  Finally,  since  the  set  DC2  is  contained  in  the  domain  of  exponential  stability,  it  follows  that  the  origin  is 
an  asymptotically  stable  equilibrium  of  the  system  (31),  with  domain  of  attraction  which  includes  the  level 
set  flC3.  In  particular,  it  can  be  shown  that  there  exists  a  function  T 2(z,  £)  which  is  continuous  and  positive 
definite  on  DCs ,  such  that 


W3  <  -T 2(z,  C)  +  2 C'PG2(M  -  Md) 


(33) 


for  all  (z,  C)  G  DCs. 


D.  Adaptive  Controller  for  the  Pitch  Moment 

The  final  step  in  our  design  is  to  determine  the  control  law  for  5e  in  such  a  way  that  | M(t)  —  Md(t)\  — >  0. 
To  this  end,  let  us  rewrite  the  thrust  and  moment  expressions  in  matrix  form.  In  particular,  let 

T  =  O^B^a)®  -  v'T(a)0 1 
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where 


Bt(ol)  =  [<a3,  a 2,  a ,  1,  0,  0,  0,  0, 0,  0,  0,  0,  0,  0,  0, 0] 

TT(<a)  =  [0, 0, 0, 0,  —  a3,  —  a2,  —  <a,  —  1, 0, 0, 0, 0, 0, 0, 0, 0] . 

Recall  the  expression  of  the  pitch  moment 

M  =  ztT  +  q  S  c  [Cm  c?  +  C^ol  +  +  Cm]  >  • 

Since  8C  has  been  chosen  as  in  (9),  by  defining 


^i5 r  i&c  ^  L 

UM  •—  nSc  ’ 


it  is  possible  to  derive  the  following  expression  for  M 

M  =  Zt{0[BT^  -  T 'T0i )  +  6>2^2^e  -  ^2^2 


where 

$2  =  [S' cCJ- ,  S c Cm,  S cCm,  S c Cm] 

b'2(x)  =  [0,0, 9,0] 

'bO)  =  [qa2,qu,0,q]. 

The  choice  of  the  update  law  for  the  parameter  estimates,  including  those  employed  for  the  velocity  controller 
earlier  in  the  Section,  is  obtained  looking  at  the  derivative  of  the  Lyapunov  function  candidate 

W4(V,z,C)=W1(V)  +  Ws(z,C) 

along  trajectories  of  the  system  defined  by  equations  (5)  and  (31)  can  be  obtained  using  equations  (7)  and 
(33)  as 

W4  <  -k1V2  +  e'1[V^1-VB1^]-T2(z,0  +  2C'PG2(M-Md) 

<  -k4V2  -  T2(z,C)  +  e'x[V  -  V  +  K  PG2{zt(6[bt<5>  -  b^i)  +  02B2Se  -  %02  -  Md] 

<  -hv2  -  T2(z,C)  +e'l[V^>1-V  +  2C 'PG2[($ZTBT  -  zT^T)0  1  -  %02  +  92B25e  -  Md] . 


Letting 

then 


0—[01i02\  ,  B3  —  [0,  B2\  ,  ^3  —  [— QztBt  +  zt'&ti  ^2]  [^3,15^3,2] 

W4  <  -hV2  -  T 2(z,  0  +e'1[V'bi-V  Bi$]  +  2(  PG2[9' B35e  -  *'30  -  Md] . 


Choosing 

Se  =  j-l*3e  +  Md] 
tf  B  3 

where  0  is  an  estimate  of  #,  it  follows  that 


O' B3Se  -  ^30  -  Md  =  {0  -o')B36e-*'3(O-0)-Md 
=  -6  B3Se  +  V'30 

=  —62B28e  +  ^35i^i  +  ^3,2^2 


where  0  =  0  —  0.  Therefore,  the  derivative  of  the  Lyapunov  function  candidate  W4  satisfies 

w4  <  -k4V2  -  T 2(z,  0  +  9AV  $1  -  V  +  2*3,1  G'2P(]  +  £,(* 3,2  -  B25e)G2P( . 
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Table  2.  Trim  conditions  considered  in  the  simulation  case  study 

Variable  Initial  trim  condition  Desired  trim  condition 

V  7846.36  ft/s  8274.96  ft/s 

h  85000  ft  96734.4  ft 

a  0.0416  rad  0.0619  rad 


Finally,  consider  the  overall  Lyapunov  function 

W(V,  z,  C,  01,  e2)  =  W4(V,  z,  0  +  0^0!  +  0'2r^02 

which  accounts  for  parameter  estimate  errors  weighted  by  positive  definite  matrices  Ti  and  IV  It  is  easy  to 
check  that  the  time  derivative  of  W  along  the  trajectories  of  the  system  satisfies 

w  =  w4  +  +  02v2 l02 

<  - kiV 2  -  T  2(z,  C)  +  0i[V  ^1  -  V  Bi§  +  24-3,1  G’2P(  +  rf§i]  +  02[(*3,2  ~  B25e)G2P(  +  Y2H2]  ■ 
By  choosing  the  update  laws  as 

0i  =  r1[F(B1$-4-1)W24'3,iG;PC] 
k  =  T2[(B25e-y3t2)G'2P(} 

one  obtains 

w  <-kiV2  -r2(z,(). 

Standard  arguments28  imply  that  all  trajectories  of  the  closed-loop  system  are  bounded,  and  that  the  tracking 
errors  V(t)  and  y(t)  are  regulated  to  zero  asymptotically. 

IV.  Simulation  Results 

To  validate  the  controller  derived  in  the  previous  section,  several  simulations  have  been  performed  on  the 
CFM  model  implemented  in  simulink®.  As  a  representative  case  study,  the  vehicle  is  initially  at  the  trim 
condition  given  in  the  first  column  of  Table  2.  The  reference  trajectory  brings  the  vehicle  to  the  new  trim 
condition  shown  in  the  second  column  of  Table  2.  The  reference  commands  have  been  generated  by  filtering 
step  increments  in  velocity  and  altitude  by  a  second-order  pre- filter  with  natural  frequency  ujf  =  0.03  rad/s 
and  damping  factor  Q  =  0.9.  The  tracking  reference  for  the  flight  path  angle  is  generated  using  equation 
(8),  where  the  gain  has  been  set  equal  to  k h  =  0.01.  The  initial  conditions  of  the  plant  parameter  estimates 
have  been  randomly  selected  within  40%  of  their  nominal  values.  Finally,  the  controller  gains  are  reported 
in  Table  3.  Figure  (1)  shows  that  the  tracking  performance  in  closed-loop  for  the  velocity,  flight-path  angle, 
and  altitude,  respectively.  It  can  be  seen  that  the  tracking  error  for  the  flight-path  angle  converges  on  a 
faster  time  scale  with  respect  to  the  other  variable,  most  noticeable  the  altitude.  The  tracking  performance 
is  quite  good,  in  spite  of  parameter  uncertainty,  and  the  presence  of  the  flexible  dynamics,  neglected  in 
the  controller  design.  The  flexible  states  remain  well-behaved,  as  shown  in  Figures  2(a)-(b).  Figures  2(c)- 
(d)  show  respectively  the  time-history  of  the  state  of  the  integral  augmentation  of  the  flight-path  angle 
dynamics  and  two  selected  parameter  estimates,  specifically  those  for  the  stability  derivative  0 io  =  SCp 
and  the  coefficient  6 15  =  SC In  particular,  while  the  estimates  0(t)  do  not  converge  to  their  true  values, 
as  the  reference  trajectory  is  not  sufficiently  rich  to  ensure  persistence  of  excitation  of  the  regressor,  the 
estimates  are  well-behaved,  and  settle  to  constant  values.  On  the  other  hand,  the  steady-state  value  of  the 
integrator  converges  to  the  trim  condition  a *  given  in  the  last  row  of  Table  2.  Finally,  Figures  2(d) -(e)  show 
respectively  thrust  and  pitch  moment,  and  the  main  control  inputs.  The  choice  of  a  relatively  low-bandwidth 
prefilter  and  a  careful  selections  of  the  controller  gains  help  maintaining  both  <F(t)  and  Se(t)  remain  within 
feasible  limits. 
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Table  3.  Design  parameters  and  gains  value 


Gain 

Value 

Gain 

Value 

kh 

0.01 

a0 

5 

h 

80 

ai 

15 

k2 

100 

Ti 

10  X  /i6xl6 

k3 

4 

r2 

10  X  /4x4 

V.  Conclusion 

In  this  paper  we  have  presented  the  design  of  a  robust  and  adaptive  nonlinear  controller  for  the  longi¬ 
tudinal  motion  of  a  air-breathing  hypersonic  vehicle.  In  an  attempt  to  design  controllers  which  are  robust 
to  model  uncertainty,  the  methodology  proposed  in  this  paper  departs  from  the  use  of  dynamic  inversion, 
in  favor  of  an  approach  based  on  sequential  loop-closure.  To  this  end,  a  variety  of  techniques  ranging  from 
adaptive  control  to  robust  semi-global  gain  assignment,  has  been  utilized  in  the  definition  of  appropriate 
intermediate  control  laws  at  each  step  of  the  design.  Simulation  results  reported  for  a  specific  case  study 
show  the  effectiveness  of  the  methodology. 

Despite  being  encouraging  in  many  ways,  the  results  presented  here  are  to  be  considered  preliminary, 
and  point  to  much  further  work  ahead.  First  of  all,  the  presence  of  flexible  effects  must  be  taken  into 
account  directly  at  the  level  of  control  design,  or  at  least  their  effect  on  closed-loop  performance  evaluated 
quantitatively.  Second,  extensive  simulations  should  be  performed  on  the  basis  of  the  original  high-fidelity 
first-principle  model,  to  validate  the  appropriateness  of  the  curve-fitted  model  for  this  type  of  design.  Since 
the  adaptive  controller  tends  to  have  a  relatively  high  dimensionality,  issues  of  controller  order  reduction 
should  be  addressed  to  obtain  more  agile  controllers.  Finally,  the  important  issues  of  limitation  in  the  control 
authority  should  be  addressed  directly  at  the  level  of  controller  design,  to  allow  more  aggressive  maneuvers, 
and  increase  the  speed  or  response  to  set-point  commands. 
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Altitude  T racking  [ft]  Flight  Path  Angle  Tracking  [rad] 


(a)  Velocity  Tracking 


(b)  Velocity  Tracking  Error 


(c)  Flight-Path  Angle  Tracking 


(d)  Flight-Path  Angle  Tracking  Error 


(e)  Altitude  Tracking 


(f)  Altitude  Tracking  Error 


Figure  1.  Simulation  Results  for  the  Control-Oriented  Model:  Controlled  Outputs. 
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(a)  Flexible  States  rn(t) 


(b)  Time-Derivative  of  the  Flexible  States,  r)i(t) 


(c)  Integrator  State  £(£) 


(d)  Estimates  of  0\o  =  SC g  and  #15  =  SC ^ 


(e)  Thrust  T(t)  (top)  and  Moment  M(t)  (bottom) 


(f)  Control  Inputs  <F(t)  (top)  and  Se(t)  (bottom) 


Figure  2.  Simulation  Results  for  the  Control-Oriented  Model:  Flexible  States,  Parameter  Estimates,  and 
Control  Inputs. 
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